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Abstract. A visualization scheme for quantum many-body wavefunctions is 
(-H ■ described, which we have termed qubism. Its main property is its recursivity: 

increasing the number of qubits reflects in an increase in the image resolution. 
Thus, the plots are typically fractal. As examples, we provide images for the 
ground states of commonly used Hamiltonians in condensed matter and cold atom 
physics, such as Heisenberg or ITF. Many features of the wavefunction, such as 
magnetization, correlations and criticality, can be visualized as properties of the 
images. In particular, factorizability can be easily spotted, and a way to estimate 
the entanglement entropy from the image is provided. 



PACS numbers: 05.30.Rt 03.65.Ud 05.45.Df 

1. Introduction 

1.1. Motivation 

Most of the difficulty of quantum many-body physics stems from the complexity of its 
fundamental mathematical objects: many-body wavefunctions and density matrices. 
In the simplest case, where we have N qubits, a wavefunction (pure state) can be 
considered as a function mapping {0, 1}^ i— > C. Therefore, it is characterized by 2^ 
complex parameters. Density matrices (mixed states) have even greater mathematical 
complexity, mapping {0, 1}^ x {0, 1}^ i-?- C, i.e. 2^^ complex parameters. 

The aim of this work is to describe a pictorial representation of quantum many- 
body wavefunctions, in which a wavefunction characterizing a chain of N qubits 
maps into an image with 2^/^ x 2^^^ pixels. Thus, an increase in the number of 
qubits reflects itself in an increase in the resolution of the image. These images 
are typically fractal, and sometimes self-similar. Extension to higher spin qudits is 
straightforward, and is also explored. Some physical properties of the wavefunction 
become visually apprehensible: magnetization (ferro or antifcrromagnetic character), 
criticality, entanglement, translation invariance, permutation invariance, etc. 
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1.2. Historical review 

Visualization of complex data is a common problem in many branches of science and 
technology. Let us review here some of the relevant hallmarks that preceded our work. 

Historically, it can be argued that the single most relevant advance in calculus 
was the discovery of the relation between algebraic functions and curves in the plane in 
the XVII century. Function visualization provided an insight which guided most of the 
subsequent development of calculus, not only by helping solve established problems, 
but also by suggesting new ones. With the advent of the new information technologies, 
complex data visualization has developed into a full-fledged field of research. The 
reader is directed to [1] for a recent review of state-of-the-art techniques, and [2] for a 
historical perspective. 

As a relevant example, the problem of visualization of DNA and protein sequences 
was addressed in 1990 by Jeffrey making use of the so-called chaos game representation 
(CGR) [3]. DNA sequences are long, highly correlated strings of four symbols, 
{A, C, G, T}. Let us label the four corners of a square with them. Now, select 
the central point of the square and proceed as follows. Pick the next symbol from 
the string. Find the point midway between the selected point and the corner which 
corresponds to the symbol. Mark that point, and make it your new selected point. 
If the sequence is genuinely random, the points will cover the square uniformly. 
Otherwise, patterns will emerge, very often with fractal structure. The original 
purpose of the technique was mere visualization, but it evolved [4] to provide 
quantitative measurements, such as Shannon entropies, which help researchers to 
characterize DNA and protein sequences iSj. 

In 2000, Hao and coworkers [B] developed a different representation technique 
for long DNA sequences that also had fractal properties. Given a certain value 
of N , they computed the frequency of every subsequence of length N within the 
global sequence, thus obtaining a mathematical object which is similar to a many- 
body wavefunction, only mapping from {A, C, G, T}^ i— ;> R. The number of different 
subsequences of length N is 4^. Hao and coworkers represented the subsequence 
probability distribution by dividing a unit square in a recursive way, into 4^ small 
squares, and attaching a color to each of them. The resulting images have fractal 
appearance, as remarked by the authors, but their quantification is not pursued. Their 
purpose is to identify which types of subsequences are under-represented, and to this 
end they analyse the corresponding patterns of low frequency. 

In 2005 Latorre [7], unaware of the work of Hao et al., developed independently 
a mapping between bitmap images and many-body wavefunctions which has a similar 
philosophy, and applied quantum information techniques in order to develop an image 
compression algorithm. Although the compression rate was not competitive with 
standard jpeg, the insight provided by the mapping was of high value [8] . A crucial 
insight for the present work was the idea that Latorre's mapping might be inverted, 
in order to obtain bitmap images out of many-body wavefunctions. 

Focusing on quantum mechanics, the simplest visualization technique is provided 
by the representation of a qubit as a point on a Bloch sphere. Early work of Ettore 
Majorana [9] proved that a permutation-invariant system of N spins-1/2 can be 
represented as a set of n points on the Bloch sphere. This Majorana representation 
has proved very useful in characterizations of entanglement [TUl [TT] . 

A different approach that can provide visualization schemes of quantum many- 
body systems was introduced by Wootters and coworkers in 2004 [12]. The idea is 
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to set a bidimensional array of operators which fulfih certain properties, and measure 
their expectation values in the given state. Those values, displayed in a 2D lattice, 
generate a discrete analogue of a Wigner function. 

1.3. Plan of this work 

In this work, we describe a set of techniques which provide graphical representations 
of many-body wavefunctions, which share many features with the schemes of Latorre 
and Hao and coworkers. The main insight is that the increase in complexity as we 
add more qubits is mapped into an increase in the resolution of the corresponding 
image. Thus, the thermodynamic limit, when the number of qubits tends to infinity, 
corresponds to the continuum limit for the images. The scheme is recursive in scales, 
and this makes the images look fractal in a natural way. In fact, as we will discuss, 
exact self-similarity of the image implies that the wavefunction is factorizable. 

In section [5] we describe the basic wavefunction plotting scheme, while section 
O is devoted to providing several examples (Heisenberg, ITF, Dicke states, product 
states, etc.) emphasizing how physical features map into plot features. The procedure 
is generalized in section 21 and some alternative plotting schemes are described, which 
allow us to try states of spin-1 systems, such as the AKLT state. Section^ on the other 
hand, deals with the fractal properties of the plots and extracts useful information from 
them. In section [5] discusses how to recognize entangled states in a wavefunction plot, 
along with a simple technique to estimate entanglement by inspection. A different 
plotting scheme, based upon the frame representation and related to the Wootters 
group approach is succintly described in section [71 and a few pictures are provided 
for the sake of comparison. The article finishes with conclusions and a description of 
future work. 

2. 2D-plot of many-body wavefunctions 

Let us consider a couple of qubits. The tensor basis is composed of four states: |00), 
|01), 1 10) and Consider also a unit square, [0, 1] x [0, 1], and divide it into four 

"level-1" squares. We can associate each of the basis states to one of the squares, as 
shown in figure [T] (top). 

The basic mapping is, therefore: 

00 Upper left 01 ^ Upper right 
10 Lower left 11 ^ Lower right 

The splitting of squares can be iterated, obtaining level-2 squares, etc., as it is 
shown in figure [T] (bottom). For a wavefunction with N qubits, we will have to descend 
down to level- squares. Each of them will correspond to one of the tensor basis 
states. If the number of qubits N is odd, the same scheme can be applied with a 
rectangular plot. The last step is straightforward: attach a color, or a gray level, to 
each of the level- iV/2 squares, depending on the value of the wavefunction. Obviously, 
using only levels of gray (or color intensity) , only real values can be attached easily to 
each tensor basis state. In order to show phases, we recourse to a color-cycle scheme. 

Figure |2] shows some features of this mapping. The ferromagnetic (FM) states, 
0000... and 1111... correspond, respectively, to the upper-left (NW) and lower-right 
(SE) corners of the image, while the Neel antifcrromagnctic (AF) states correspond 
to the other two corners: 0101... is the upper-right (NE) corner, and 1010... is the 
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Figure 1. 2D plotting scheme of many-body wavefunctions. Top: each of the 
tensor basis states for 2 qubits; |00), |01), |10) and |11) is mapped into one of the 
four level-1 squares. Bottom: mapping of 4-qubit basis states into level-2 squares. 
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Figure 2. Exemplification of some features of the 2D-pIot scheme. The FM 
states 0000... and 1111... correspond to the NW and SE corners of the image, 
respectively. The NE and SW corners, on the other hand, correspond to the Neel 
AF states. The Z2 symmetry operation 1 corresponds to a rotation by 180° . 



lower- left one (SW). It is straightforward to realize that the Z2 symmetry operation 
-H- 1 corresponds to a rotation of 180° around the center of the plot. 

Let us consider any state s E {0, 1}^ and denote its bits by s = 
{X1Y1X2Y2 ■ ■ ■ XnYn}, with n = N/2. In order to find the point in the unit square 
where this state will be mapped, build the following numbers: 



n n 



(2) 
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(A) (B) (C) (D) 

Figure 3. Qubistic plots of some important states with N = 4 qubits. (A) The 
simplest factorizable state |0000). (B) The GHZ state. (C) W state. (D) Dicke 
state at half-filling. 

Those are the coordinates of the upper- left corner of the corresponding level-n square, 
if (0, 0) is the upper-left corner of the square, and the y-coordinate grows downwards. 

In our plots, unless otherwise stated, each cell is filled with a color corresponding 
with its wave-function amplitude according to the following scheme: color intensity 
corresponds to the modulus (white means zero), and hue is used as a phase indicator. 
Concretely, red is used for positive values and green for negative ones, with a smooth 
interpolation scheme. 

Figure [3] provides some simple examples of states with = 4 qubits. Panel (A) is 
the qubistic plot for the factorizable state |0000), in which only the upper- left corner 
cell of the plot is marked. In (B) the GreenbergerHorneZeilinger (GHZ) state is shown, 
|0000) -I- In this case, two opposite corner cells are marked, with the same color 
since their relative phase is positive. The third panel, fig. |3] (C) corresponds to the 
so-called W state for TV = 4 qubits, i.e.: |1000) -I- |0100) -I- |0010) + |0001). The plot 
consists of N marked cells distributed along the upper and leftmost rows of the plot. 
For larger values of N the spacing among these marked cells becomes exponential. The 
fourth plot in figure [3] is the Dicke state at half-filling, i.e.: the linear combination, 
with equal weights, of all basis states with half the qubits 1. 

3. Examples of Qubistic 2D-plots 

Along this section we will study qubistic plots of low-energy states of hamiltonians 
which are relevant in condensed matter physics and ultracold atomic cases, giving 
special attention to quantum phase transitions (QPT) [T51 [T^ . 

3.1. Heisenberg Ground State: Spin Liquid structure 

Our next example will be taken from the low-energy spectrum of the antiferromagnetic 
(AF) spin-1/2 Heisenberg model in ID with periodic boundary conditions (PBC). 

N 
i=l 

The top panel of figure 0] shows the ground state of equation [31 while the bottom 
row shows the first three excited states, which constitute a spin-1 triplet. 

Let us focus on the ground state (figure 01 top panel). The most salient 
feature is its intense diagonal line, joining the two Neel states, which get maximal 
weight. The states conforming that diagonal are all made up of pairs 01 and 10, 
in any order. This main diagonal is the depiction of a set of pairwise singlet bonds: 
(1,2)(3,4)...(A-1,A). 
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Figure 4. Top: Ground state of the ID spin-1/2 AF Heisenberg model with 
PBC for 12 spins; Bottom: the lowest energy excitations, which make up a 
triplet. White means zero probability, color intensity reflects the modulus of 
the wavefunction amplitude, while color hue marks the phase. 



There is another interesting feature in this image. The two parahel diagonal hnes 
with slope 1/2 have the same intensity as the main diagonal. What is their origin? 
A clue can be obtained when we depict the GS of the Heisenberg model with open 
boundary conditions (see figure [5]) . It is apparent that these secondary lines have 
almost disappeared. In order to finally clarify the nature of these secondary lines, let 
us consider 7?., the right-shift translation operator (with periodic boundary conditions). 
If TZ acts on the states composing the main diagonal, the result is the two secondary 
diagonals, and viceversa, as it can be seen in figure HI It is now straightforward to 
provide a physical interpretation: the secondary diagonals depict the other possible 
set of pairwise singlet bonds: (2, 3)(4, 5) • • • {N, 1). When periodic boundary conditions 
are employed, both structures are equally important, but not under open ones. 

The slope 1/2 of those secondary diagonals can be understood as follows. 
According to equation [2J acting with the right-shift translation operator 7?, on a state 
given by bits {X1Y1X2Y2 ■ ■ ■ X„Yn} we obtain {y„XiYiX2 • • • r„_iX„}. Thus, 7^ 
maps the point {x,y) into a point very close to ((y -|- Yn)/2,x). Consequently, the 
image of the x — y line is approximately x = {y -\- K„)/2, i.e.: the two secondary lines. 
A second application of the right-shift operator TZ on these two secondary lines returns 
the original main diagonal. Of course, the same effect is obtained with a left-shift. 
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Figure 5. Ground state of the ID spin-1/2 AF Heisenberg model with open 
boundary conditions and N = 12 spins. Notice that, as opposed to the case of 
figure |4] the secondary diagonals have almost vanished. 




n 



Figure 6. The two diagonal structures which make up the Heisenberg ground 
state are related through a right-shift translation operator TZ.. If denotes a 

singlet bond, on the left we have (1, 2){3, 4) ■ ■ ■ (A^-l, A^), and (2, 3)(4, 5) ■ ■ ■ {N, 1) 
on the right. 



3.2. Next-nearest-neighbour Heisenberg: Marshall rule and Frustration 

Still there is one more interesting feature of the image of the ground state of the 
Heisenberg Hamiltonian, |31 According to Marshall's rule [H], the sign of each 
wavefunction component of the ground state of the Heisenberg AF model in a bipartite 
lattice (split into sublattices A and B) can be given as (—1)^-* , where Na is the number 
of up-spins in sublattice A. In our case, a ID lattice with PBC, the two sublattices 
are just the odd and even sites. It is not hard to recognize that, if we select the odd 
sites to make up sublattice A, then the sign rule tells us that all states in the same 
horizontal line must have the same sign. But, on the other hand, if sublattice A is 
made up with the even sites, then the rule tells us that all states in the same vertical 
line will have the same sign. Both conditions can be fulfilled, both in the PBC and 
the OBC figures, H and m 

Marshall sign rule can not be applied if the system presents frustration, i.e.: 
when the Hamiltonian couples spins in the same sublattice A B. Let us consider 
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Figure 7. Destruction of the sign pattern given by Marshall rule as J2 increases 
for the ground state of the Hamiltonian given in eq. U, with N = 12 qubits 
and PBC. From top to bottom, the values of J2 are J2 = 0, 0.2 and 0.3. The 
left column shows the full image. The right one only depicts the phases of the 
non-zero amplitude values. 



the next-nearest-neighbour AF Heisenberg Hamiltonian (also known as J1J2 model): 

AT AT 

H=JlJ2S^■S.+l + J2Y,S^■S,+2 (4) 

1=1 1=1 

where Ji = 1 and J2 > 0- Then, as J2 increases, the system undergoes a quantum 
phase transition (QPT) at around J2 ~ 0.24. Figure [7] shows how the sign-structure 
is destroyed slowly as J2 is increased from J2 = to J2 = 1/2- The point J2 = 1/2 
is special, since the ground state is then exactly known: the Majumdar- Ghosh point. 
Its rather simple structure is apparent in figure HI 
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Figure 8. Majumdar-Ghosh state, ground state of the Hamiltonian in equation 
Hwith Af = 12 qubits and PBC for J2 = 0.5. 

3.3. Product States 

Let us now consider the simplest possible quantum many-body wavefunction: a 
translationally invariant product state, defined as 

|vl/) = (a|0) + /3|l)f (5) 

Factorizability is a very strong property, which shows itself in a very appealing 
way in our plots. Figure IHl shows such a product state in the cTz basis. Physically, 
factorizability implies that measurements performed on any qubit should have no 
influence on the remaining ones. Concretely, we can measure az on the first two 
qubits. If the result is 00, the wavefunction-plot which describes the rest of the 
system will be (a normalized and rescaled version of) the upper-left quadrant of the 
plot. Correspondingly, if the results are 01, 10 or 11, the wavefunction-plot will be just 
a (normalized and rescaled version of) other quadrant. Thus, factorizability implies 
that all four quadrants are equal (modulo a normalization factor). This line of thought 
can be extended to the set of the first 2k qubits, thus showing that if we split the plot 
into a 2*^ X 2*^ array of sub-images, they should all coincide, modulo a normalization 
factor. This gives the characteristic look to the plots of product states. We will return 
to this topic in section [HI when we discuss entanglement. 

3.4. Dicke states 

Another interesting example is provided by the so-called Dicke states [16 . Those are 
defined as the linear combination, with equal weights, of all tensor basis states with the 
same number rig of I's in their decomposition. In our examples, we will focus on the 
half-filling case, rig = N/2. In fermionic language, they constitute the ground state of 
a free-fermion model with homogeneous difl^usion on a complete graph at half-filling, 
and in spin-1/2 language it is the Sz = Q component of the maximal spin multiplet. 
Figure [10] shows the pattern obtained for different lattice sizes. It is apparent how 
a fractal develops. Their similarity to the right column of figure [7] is, of course, not 
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Figure 9. Product state following equation |5] with /3/a = 1.3. 



■ 

J 










r 








. .J 




















r 


fr r ■ 




Figure 10. Half-filling Dicke states for N - 
the fractal structure develops. 



8, 10, 12 and 14 qubits. Notice how 



casual: the ground states of the Heisenberg-like models have global magnetization 
zero, which make them similar to half-filling Dicke states. 
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3.5. Ising model in a transverse field: Criticality 

As a different relevant example, let us consider the spin- 1/2 AF Ising model in a 
transverse field (ITF), in a ID chain with PBC: 

N N 
i=l 1=1 

For F = 1, the system presents a quantum phase transition (QPT). Figure [TT] 
shows the plots obtained from the GS for different values of F. For F — 0, the ground 
state consists only of the two Neel states. As F increases (first two top panels), the 
points which come up first correspond to a single defect, at all possible positions in the 
lattice. The non-zero probability amplitudes extend further away from the original 
corner states as F approaches criticality, and at that point Fc = 1, the non-zero values 
have extended through the whole image, albeit quite inhomogeneously. From that 
point, increasing F makes the image more and more homogeneous. For F — > oo, the 
ground state would consist of all spins pointing in the A"-direction, and this implies 
that the wavefunction components will all take the same value. 

3.6. Infinite-range Hamiltonians 

And let us finish this section by considering infinite-range Hamiltonians, i.e.: those in 
which all spins are linked to all others. They can be thought of as infinite-dimension 
or mean-field systems. Those Hamiltonians commute with the full set of generators 
of the permutation group. Therefore, their ground states are often invariant under 
it. Compared to translation invariance, this symmetry group is so large (TV! elements 
vs. N) that it leaves very little freedom: a fully permutation-invariant wavefunction 
of N qubits is characterized by just TV -I- 1 independent components, one per global 
magnetization sector. Thus, permutation-invariant wavefunctions have a clear visual 
fingerprint. 

Figure[T2 (left) shows the GS of the infinite-range AF ITF Hamiltonian for F = 1 
and iV = 12 qubits, illustrating this high degree of symmetry. The right part of the 
figure shows a random permutation-invariant state. It is not a coincidence that it 
reminds so strongly of the Dicke states, since each magnetization sector shares the 
same color. The infinite-range Heisenberg Hamiltonian ground state is not shown 
because it is strongly degenerate, so invariance under the permutation group remains 
only as a property of the full subspace. 

4. Other Plotting Schemes 

4.1. General Formulation 

The previous procedure can be generalized in the following way. Let D be any 
domain in R*^, which can be partitioned into m congruent subdomains SiD, with 
i G {0, • • • ,TO— 1}, all of them similar to D. In our current example, D — [0, 1] x [0, 1], 
the unit square, which is divided into m — A smaller squares, which we denote by SqD 
(upper-left), SiD (upper-right), S2D (lower-left) and S^D (lower- right). 

The action of operators Si can be iterated. Thus, S1S3D denotes the upper- 
right quadrant of the lower-right quadrant of the original square. We can define a 
geometrical index as a sequence of integers Ig = {*fc}fc=i: with ik G {0 • • • m — 1}. 
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Figure 11. Ground states of the AF ITF Hamiltonian with A'' = 12 qubits 
and PBC. Values of the transverse field are: (top) P = 0.2 and 0.6, (central) 
r = = 1, (bottom) F = 1.4 and 1.8. 

Each geometrical index denotes a (small) domain Si-^ ■ ■ ■ Si^D, similar to the original 
one. In our example, a tiny square. We can, thus, define a mapping S which converts 
geometrical indices into regions of M'^ which are similar to D: S{Ig) = Si^ ■ ■ ■ Si^D. 

Now let us focus on the tensor-product structure of the quantum Hilbert space. 
Each state is characterized by a quantum index, i.e.: a set of indices taken from 
certain discrete finite set: Iq G E^. In our case, E ~ {0, 1}. In the case of spin-1 
systems, E = {—1,0,1} or, more simply, E = { — ,0,+}. 

The last piece of the scheme is a function A4 mapping quantum into geometrical 
indices, M : E^ i-^ {0, • • • , m— 1}", such that Iq = M{Iq). In our case, this function 
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Figure 12. Left: Ground state of the AF infinite-range ITF model for F = 1 witii 
JV = 12 qubits. Riglit: a random permutation-invariant state, also with A'' = 12 
qubits. 

groups the quantum indices in pairs, and combines each pair into a single geometrical 
index with the simple binary mapping: 00 — 0, 01 — )■ 1, 10 — >■ 2, 11 — >■ 3. It should be 
noted that n = N/2. This mapping should be bijective, so as not to lose information. 

Now, the full wavefunction plotting scheme /C is defined by providing the original 
region, D, the set of similarity transformations, {Si}, i G {0, •••,m — 1} and the 
indices mapping function A4. Thus, /C(/q) will denote the region in R'' obtained by 
applying S to the geometrical index associated to Iq, i.e.: /C(/q) = iS(A^(/q)). Those 
cells make up a partition of D. It is easy to prove the essential properties: 

Uq6E"/C(/q) = D 

/c(/Q)n/c(/^)-0 ^ Iq^i'q (7) 

Thus, for every x G £>, there exists a single Iq e such that x G JC{Iq). 
This property ensures that we have can pull-back wavefunctions, i.e.: functions 
ip : I—)- C, into complex-valued functions on D, K.{ip) : D i-^ C 

So, can we devise other possible plotting schemes? Will they make different 
properties apparent? We will approach those questions in the rest of this section. 

4.2. ID plot 

The simplest possible plotting scheme can be realized in ID for qubits. Let D be the 
[0, 1] segment, split every time into two halves: So selects the left part, and the right 
one. Now, the resulting /C mapping is equivalent to a binary lexicographical ordering 
of the wavefunction components. More explicitly: divide the domain [0, 1] into 2" 
equal cells, index them from to 2" — 1 and attach to each of them the wavefunction 
component with the same associated index. 

Figure [13] shows plots (ID) of the ground state of the antiferromagnetic ITF 
model, equation [6l for several values of F. This plotting scheme is, evidently, much 
less appealing than the bidimensional ones. On the other hand, its simplicity is helpful 
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when attempting to clarify some of the features, e.g. the Fourier analysis made in 
section [5l 

4.3. Spin- 1 plots: AKLT states 

As we have already stated, the quantum indices can be built upon any local set of 
quantum numbers. For a set of spin-1 particles, the choice is E = {—1, 0, 1} or, more 
simply, S = { — ,0,+}. If we start with the same domain, D — [0,1]^, the natural 
decomposition is into 3x3 subdomains, as shown in the following (cartesian product) 
scheme: 





-0 


++ 


0- 


00 


0+ 


+- 


+0 


++ 



(8) 



Of course, this is not the only possible mapping. With this one, we have chosen 
to show the structure of the Afheck-Kenedy-Lieb-Tasaki (AKLT) state [1^ . It is the 
ground state of the following Hamiltonian: 

N 1 

H = Y.^^- '^'■+1 S^+lf (9) 

1=1 

This state is an example of valence bond solid (VBS), and has attracted 
considerable attention because of its relation to the Haldane conjecture (TB], its non- 
local order parameter [TH] and as a source of inspiration of tensor- network states [50] . 

The result can be seen in figure [TJ] where, for better visualization, we have marked 
only the non-zero components of the wavefunction. Notice the strong self-similarity 
appearance. 
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Figure 14. Ground state of the AKLT spin-1 Hamiltonian, for TV = 6, 8 (top) 
and 10 (bottom) spins. Notice how the fractal structure develops. 

4.4- Alternative Square Plot 

Restricting ourselves to qubits and V ~ [0, l]'^, it is still possible to have another 
inequivalent plotting scheme, by changing the assignments: 

00 Upper left 01 Upper right , , 

11 Lower left 10 ^ Lower right 

In this new plotting scheme the two left corners (top and bottom) represent the 
FM states, and the right corners the Neel states. 

It can be shown that these two are the only possible inequivalent plotting schemes 
for qubits on [0, 1]^. The reason is the following. There are 4! = 24 possible 
associations between {00, 01, 10, 11} and the four quadrants. The group of symmetries 
contains three rotations, two reflections on the horizontal and vertical axes and two 
reflections on the two diagonals, i.e.: 12 different elements. This leaves only 4!/12 = 2 
inequivalent choices. 

As an example, figure [15] (left) shows the ground state of the critical (F — 1) 
ITF model with TV = 12 qubits (eq. [6]). It is therefore, an alternative pictorial 
representation of figure[TlJ Figure [T5] (right) shows the ground state of the Heisenberg 
model with iV = 12 qubits, in the new plotting scheme. The Neel states are 




Figure 15. Alternative 2D representations, following equation 1 101 with A'^ = 12 
qubits and PBC. Left: Ground state of the critical AF ITF model. Right: Ground 
state of the Heisenberg model. Compare, respectively, to figures [TT] (central) and 

HI 



now situated in the lower and upper right corners. Therefore, the main diagonal 
line, hallmark of the spin-liquid structure, lies now in the rightmost vertical line. 
The secondary diagonals, on the other hand, are now dispersed, in a Sierpiriski-like 
structure. 

It is apparent that figure [15] (bottom) is smoother than its counterpart, figure [TT] 
(central). All possible plotting schemes are equally valid, in principle, just as a polar 
and a cartesian representation of the same function are. Can we provide some sense of 
plotting quality'? Perhaps: a smoother plot suggests that the neighbourhood structure 
of the original wavefunction is respected more properly by the plotting scheme. 



4-5. Triangular Scheme 

It is possible to design a 2D plotting scheme of qubits which does not require grouping 
the quantum indices in pairs. Let Dhe a, rectangular isosceles triangle of unit side, with 
vertices at (—1,0), (0, 1) and (1,0). It can be split into two similar triangles, of side 
l/-\/2. Let So and be the operators which select the left and right triangles (as seen 
when the right-angle vertex is up). Figure [16] shows how such a representation maps 
bits into cells. Within this scheme the Neel states go towards the left and right bottom 
corners. The FM states correspond to two points near the center, symmetrically placed 
with respect to the height. In the thermodynamical limit, these FM points can be 
obtained summing a geometrical series: (±1/5,2/5). 

Figure [TT] depicts the ground states of the critical ITF and the Heisenberg model, 
and a product state. Notice that the diagonal lines in the original representation 
for the Heisenberg GS have mapped now to the perimeter of the triangle. The main 
diagonal is the hypotenuse, the two secondary diagonals are the other two sides. The 
remaining structure, which was not quite clean in the original representation, comes 
here as a Sierpihski-like structure. 
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Figure 16. Plotting scheme using rectangular triangles. Top and center; first 
two iterations. Bottom: First four iterations (dashed lines). The two Neel states 
of the fifth iteration are the bottom left and right filled corner triangles. The two 
FM states correspond to the symmetrically placed filled triangles near the center. 



4- 6. More Exotic Plotting Schemes 

We will now propose other plotting schemes in order to show the versatility of the 
procedure. 

In the case of spin-1 systems, the only alternative that we have found in order 
to make the quantum and the geometrical indices coincide is to work on a Sierpinski 
triangle. The original domain is, this way, naturally split into three similar domains: 

5- D, SoD and S+D. Nonetheless, it has the disadvantage that the domain is not 
simply connected. 

Even more exotic plotting schemes are conceivable. Let Aq be a regular hexagon. 
Now proceed to build Ai as the union of a-nd six congruent hexagons built upon its 
sides. Repeating the scheme, and rescaling at each step, we reach a fixed point: Aqo, 
with the following property: it can be split naturally into 7 similar cells of exactly the 
same shape [H]. 



]ubism: self-similar visualization of many-body wavef unctions 



18 























y 


y 










y 




y 






y y 


^ y 


y k> 
> y 


> y 


> y > 



Figure 17. Triangular representations of many-body wavefunctions. Top: GS of 
the AF-ITF model with = 12 qubits and PBC. Center: GS of the Heisenberg 
model with N = 12 qubits and PBC. Bottom: product state with N = 12. 



5. Self-similarity of the wavefunction plots 

The plotting schemes described in the present article are evidently self-similar. It 
is obvious that the first qubit determines the largest-scale properties of the plot, and 
subsequent qubits determine lower scales properties. The question that we will address 
is: how does this self-similarity of the scheme map into fractal or self-similar properties 
of the wavefunction plots? 

5.1. Translation-invariance and self- similarity 

Let TZ be the cyclic right-translation shift. A wavefunction has translational symmetry 
if, for any quantum index /q, |-0(/q)| = |'0(7?.(/q))|. Does this symmetry bear any 
visual consequences in the wavefunction plots? 
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Figure 18. Illustrating translation invariance properties in a qubistic plot. All 
three plots represent the qubistic representation of a Af = 6 wavefunction. Assume 
we measure cr^ on qubits 1 and 2, and the results are and 0. Then, after 
the measurement, the wavefunction-plot will be given by the marked cells of the 
top plot. The central and bottom plots are equivalent, but with measurements 
performed on qubits {3, 4} and {5, 6} respectively. If the wavefunction is 
translation-invariant, all three resulting wavefunction-plots should be exactly 
equal. 

In a translationally invariant system, a measurement performed on the first two 
qubits and another performed on the last two should have the same effects. Let us focus 
on a given possible outcome of the measurement, e.g.: 00. Now, the wavefunctions 
describing the rest of the system should coincide. If the measured qubits are the first 
two, the new wavefunction-plot is obtained by selecting the upper-left quadrant of the 
original plot. On the other hand, if the measurement has been done on the last two 
qubits, we should decimate: group the plot pixels into 2x2 blocks, and select the 
upper-left pixel out of each block. Both wavefunctions should coincide, as a result of 
translation- invariance. 

So, plots of translation-invariant wavefunctions display self-similarity in the 
following sense. Divide the plot into a matrix of 2^ x 2^ sub-plots (A; e {0, • • • ti ~ 1}) 
and do a further division of each sub-plot into 2x2 quadrants. Selecting the same 
quadrant from each sub-plot and rebuilding a full image will yield the same result, for 
all possible values of k. Figure [TSl illustrates the criterion. 
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Figure 19. Fourier transform of the ground state of the ID plot of the critical 
ITF model, with Af = 10 qubits and PBC, for T = 0.6. 



5.2. Measures of scale invariance 

Scaling invariance of the wavefunction plots should also be visible in the Fourier 
transform. In effect, figure [19] shows the transform of a ID-plot of an AF ITF 
Hamiltonian (eq. ^ with iV = 10 qubits and PBC. The momenta are displayed 
in logarithmic scale, and we can spot a clear periodic structure. Evidently, exact 
log-periodicity is impossible to achieve since each period contains a larger number of 
degrees of freedom than the preceding one. This feature is visible for a wide range of 
transverse fields, i.e.: it is not linked to criticality. 

Another interesting indicator of self-similarity is provided by the Renyi fractal 
dimensions [?]. Let us consider the probability distribution associated to a 
wavefunction plot (taking the modulus squared), Pn — {pN,i} for N qubits. We 
can compute the Renyi entropy of order q, i.e.: 



where 6 is 2 for qubits or 3 for spin- 1. With this notation, do, di and c?2 are, 
respectively, the support, information and correlation dimensions of the fractal. The 
full set of dq provide the same information as the multifractal spectrum. 

Figure [5D] shows a few Renyi dimensions for an AF ITF model as a function of F. 
In our case, the support dimension do is always 2, since all probability values are non- 
zero. All the other dimensions interpolate between (for F — >■ 0) and 2 (F — >■ oo). The 
information dimension, di seems to capture most accurately the physical properties 
of the model, since its growth rate is maximal at the critical point. 

It is an interesting exercise to prove that, for the AKLT state shown in figure [Til 
all Renyi dimensions with g > are equal to log(4)/ log(3) sa 1.26. 




(11) 



Now we define the Renyi dimensions by: 



d„ . lim 




(12) 
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Figure 20. Renyi fractal dimensions (equation I12|l of the wavefunction plot of 
the ground state of the AF ITF model with PBC. Notice that the X-axis is labeled 
with arctan(r) in order to highlight the symmetry. The vertical bar shows the 
critical case, = 1. Computations are carried out for = 10 qubits, but the 
results are quite independent of the system size. 

6. Visualization of entanglement 

One of the most intriguing features of quantum many-body systems is entanglement. A 
system is entangled if measurement on one of its parts affects the results of subsequent 
measurements on others, even if they are well separated. Einstein himself described 
this phenomenon as "spukhafte Fernwirkung" (spooky action at a distance) [22]. It is 
considered as a resource for quantum computation and communication [23j , as well as 
providing very useful insight regarding quantum phase transitions |24) . 

6.1. Visual estimate of entanglement 

Is entanglement visualizable from our wavefunction plots? Yes. Summarizing the 
results of this section we may say that entanglement shows as image complexity. Let 
us consider all quadrants of level-fc within the plot, normalized. If there are only 
p different quadrants, then the entanglement entropy is < log(p). Concretely, if all 
level-fc quadrants are equal, the system is factorizable. 

In section [3.31 we discussed product states, i.e.: systems without entanglement. 
Let us recall the conclusions exposed in that section. If the first two qubits are 
disentangled from the rest of the system, measurements made upon them should not 
have influence on the rest. Therefore: all four quadrants of the plot are equal (modulo 
normalization). If all the qubits are disentangled (at least by pairs), then the result 
is extended: if the plot is split into a 2^^ x 2^ matrix of subimages, for all fc, all the 
sub-images are equal (modulo normalization). This result can be expressed in a more 
concise way: the plot of a product state is trivially self-similar. Every quadrant, of 
every size, is the same as any other, after proper normalization. 

What happens if the system is entangled? Let us now consider a generic 
wavefunction, j^E'), and split the system into a left and a right parts, L and R. The 
left part will correspond to qubits 1 to 2fc and the right part to qubits 2fc + 1 to iV, 
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for any k. We can always perform a Schmidt decomposition: 

rn 

4=1 

where the orthonormal sets {jV'i^)} and {jV^f^)} are cahed the left and right-states, 
and characterize the physics of each part, Xi are called the Schmidt coefficients and 
m is the Schmidt rank, which is a measure of entanglement. If rn = 1, the state is 
factorizable. A state with Schmidt rank m can not have entanglement entropy larger 
than log(m). 

The left part corresponds to the larger scales, and the right part to the smaller 
ones. Let us make this statement concrete. 

Consider the Hilbert space for the left part, and let {\x)} be the basis of tensor 
states for it. E.g.: if A; = 1, the left part has two qubits and the states {|a;)} are |00), 
|01), 1 10) and |11). Now we will consider what are their geometric counterparts in the 
wavefunction plot. Within the original 2D plotting scheme, qubits 1 to 2k correspond 
to the first k quadrant divisions. Let us divide the original plotting square into a 
matrix of 2^ x 2'' quadrants. Each tensor state \x) can be attached to one of these 
quadrants, which we will denote simply by x. 

The left-states, ) can be expressed as 

X 

Now, let us focus on the right part. Each right-state, \ipf') can be plotted inside 
a level-fc quadrant using the standard representation. Let us call the corresponding 
plot Rj. 

What is the actual image, for the full wavefunction plot, on the x-th quadrant? 
Inserting equation [HI into the Schmidt decomposition [T51 we can see that it is given 
by the expression 

m 

C{x) = ^^^-^^ (15) 

1=1 

The conclusion is that, for each level-k quadrant the plot is a linear combination 
of the m right-state plots, with weights given by the the left-states components and 
the Schmidt coefficients. 

Therefore, level-fc quadrants within the final plot are built upon only m 
fundamental images, or building bricks, which are the plots of the right-states. In other 
terms: the Schmidt rank m for a given left-right partition coincides with the effective 
dimension of the subspace spanned by all images in level-k quadrants. This statement 
provides a way to give a coarse estimate the entanglement of the wavefunction: if, at 
level-fc, the number of different quadrants is p, then the block of the first 2fc qubits 
has a Schmidt rank oi m < p, and the entanglement entropy is S < log(p). As a 
corollary, if all quadrants are exactly the same, then m = 1 and S — 0, the system is 
factorizable, as we already stated. 

The logic behind the estimate is to find the number of different building blocks 
at each scale. If we want to be precise, the Schmidt rank is given by the dimension of 
the subspace spanned by all quadrant images at a certain level, but this value is much 
more difficult to estimate visually. 
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Figure 21. Visual estimates of entanglement using qubistic plots for some 
relevant states of A'^ = 4 qubits. The first row shows the usual qubistic plots 
in the basis. The second row, on the other hand, depicts the plots using the 
CT^ basis. Since entanglement is invariant under local changes of basis, the visual 
estimate of entanglement should not change. The third row provides the Schmidt 
rank in all cases for the separation between the first and second pairs of qubits, 
and the fourth row contains the von Neumann entropy (logarithms are always 
with basis 2). The states correspond to the columns: A.- |0000); B.- GHZ state; 
C- W state; D.- Dicke state at half filling; E.- |0000) + |1111) - |1010) - |0101>. 



Let us apply the estimate in a set of simple cases, with N — A qubits (i.e.: only two 
levels). Figure [21] shows the qubistic plots for a set of states similar to those of figure[3] 
Since entanglement is invariant under local changes of basis, we also show the qubistic 
plot in the basis of eigenstates of . Both plots provide a similar estimate, which 
is compared in each case with the exact value. (A) The state |0000) is factorizable, 
which can be seen in both plots. In the tr^ plot, only one of the sub-images is non-zero. 
In the picture, all four sub-images are the same, modulo a sign. (B) The GHZ is 
not factorizable. In both basis can be seen that the number of different sub-images 
is 2. (C) corresponds to the W state, which is a bit more complex. In the ti^ basis 
it is evident that the number of different sub-images is 2, which corresponds to the 
Schmidt rank. In the cr^ basis, the visual estimate gives 3 different sub-images. Our 
prediction is still valid, since the estimate only provides an upper bound. The reason 
for the error is that the 3 sub-images are not linearly independent. This example 
serve as a warning: some basis may provide clearer visual estimates than others. 
(D) Is the Dicke state at half-filling. In this case the visual estimate coincides for 
both basis, 3 different sub-images. But do not have the same weighty and the von 
Neumann entropy is smaller than log(3). (E) The |0000) + |1111) - |1010) - |0101) 
state has four different sub-images in both basis, and achieves maximal Schmidt rank 
and entanglement entropy. 

The strategy can be applied to the AKLT state, depicted in figure [TH At any 
splitting level the exact number of different images is always 5. But, as the number 
of sites increases, some of these images become more and more alike, until only 3 of 
them are distinguishable. Figure [22] shows the sub-image pattern more clearly. See, 
for example, the — h and 0+ quadrants of the plots in figure [14] their differences are 
easy to spot for = 6, but almost unnoticed for N = 10. This implies that the 
Schmidt rank is always < 3, providing the estimate S < log(3), independent of the 
depth level, which is exactly the actual value in the thermodynamic limit [25] . 
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Figure 22. AKLT qubistic plot for N = 8 spins and tlie associated sub-image 
pattern: boxes with the same letter contain (very approximately) the same sub- 
image. 
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Figure 23. Entanglement entropy of a block of I sites in a half-filling Dicke state 
of L = 10, 12 and 14 sites, compared to the upper bound obtained from the 
qubistic plot, which is S{1) < \og{l + 1). 



On the other hand, taking the half-filUng Dicke states of figure ITOl it is evident 
that, at every magnification level, the number of different subimages increases by two. 
Thus, S{k) < log(2fc + f ) in terms of levels, or S{1) < log{l + 1) for qubits, if I < N/2. 
This bound is found to be fulfilled by the numerical calculations shown in figure [23l 

The reason for the difference between the estimate and the actual values of 
entanglement in fig. [53]is twofold. First, the number of different level-^ quadrants is, in 
general terms, a very poor way to estimate the dimension of the subspace spanned by 
them. Second, the value estimated this way is just the Schmidt rank, whose logarithm 
is just an upper bound to the actual entanglement entropy. Both problems can be 
handled within the more comprehensive framework, described in the next section. 

6.2. Entanglement and the cross-correlation matrix 

Given a wavef unction plot and a level fc, let us divide the full region into a grid of 
2'' X 2'^ sub-plot. Moreover, let x be an index running through all such sub-plots 
and C{x) be the actual image displayed in it, as in equation 1151 Now we define a 
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cross-correlation matrix for the plot image, R{x,x'), as 

R{x,x') = {C{x)\C{x')) (16) 

This cross-correlation matrix bears full information about entanglement of the 
first 2k qubits within the wavefunction, as we proceed now to show. 

According to equation 1151 the image on quadrant x is given by a linear 
combination of the right-states. Using the orthogonality property assumed for them 
we obtain 



m 

R{x,x') = ^ (V-f I V'f^A.A.Vi;' l^f ) 

J2xfi^^^g, (17) 

1=1 

Thus, we recognize that R{x, x') is just the density matrix for the left part. In other 
terms: 

Rix,x') = p^,, (18) 

Therefore, the cross-correlation matrix of the wavefunction plot holds full 
information related to entanglement. 

For example, for a product state, all subimages are equivalent modulo 
normalization. Thus, we can assume that 

Rf{x,x') ^ N{x) ■ N{x') (19) 

with N{x) — {C{x)\C{x))^/^ is the norm for each subimage. Obviously, ^^i^) — 1 
and, thus, the matrix Rp is just a projector on a line. Its spectrum is, in decreasing 
order, (j{Rf) = {1, 0, • • • , 0}. Therefore, its entanglement entropy is zero. 



7. Frame representation 

In this last section we describe a rather different approach to the problem of providing 
a graphical representation of a quantum many-body system, but still self-similar by 
design. Instead of plotting wavefunction amplitudes, or probabilities, we can plot the 
expectation values of a bidimensional array of operators, chosen in such a way that 
the full information contained in the wavefunction is preserved. This is called a frame 
representation of the quantum state |26) . According to Wootters and coworkers [1^, 
the final representation may correspond to a discrete analogue of a Wigner function 
|27j . with very interesting properties in order to characterize non-classicality, such as 
its negativity [28 . 

Let us consider a system of n qubits, described by a certain density matrix p. 
Now, let us consider the unit square [0, 1] x [0, 1] and any two numbers, x and j/, 
characterized by their binary expansion: x — Q.X1X2 ■ ■ ■ Xn, y = 0.1112 The 
value attached to the point (x, y) in the plot will be given by the expectation value in 
p of the operator A{x,y): 

f{x,y)^iY[pA{x,y)] (20) 

where A{x, y) is given by: 
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Figure 24. Illustrating the frame representation of eq. 1221 Top; operator 
assignment for a single qubit. Bottom: For two qubits. Products must be 
understood as tensor products, with the superscript denoting the qubit index. 



Aix. y) - (g)(-*)^'=^^ {alf^ {alf" (21) 

k=l 
n 

k=l 

In other words, we plot the expected value of every combination of tensor products 
of {cr" = /, (7^ = (T^, = , = a^}. In particular, on the y = line we get 
uniquely correlations in CTj,; on x = 0, those in <7y, and on x = y those corresponding 
to az- Such representation is unique for every density matrix, and can be reverted as 
follows: 

^ = E^/(^'y)^(^'y)- (23) 

x,y 

In order to attain some intuition about the representation, figure [24] illustrates 
it for one and two qubits. At each cell, we depict the expected value of a "string" 
operator, as shown. 

Figure [55] shows our first example: the frame representation of a product state 
given by 

|vl/) = (cos(j)|0)~sin(^)|l))' (24) 

i.e.: a spin pointing half-way between the —X and Z axes. The plot shows a striking 
Sierpihski-like structure, which can be fully understood by noticing that, in this state, 
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Figure 25. Frame representation for a particular product state of = 6 qubits, 
described in equation 1241 Notice the Sierpinski-Iike structure, whicii is explained 
in the text. 



{^\ I*) and (*! 0-2 I*) are nonzero, while {^\ Oy = 0. If, in figure [M] (bottom) 
we cross out all elements with a cr^, the Sierpihski-like structure will appear. Self- 
similarity, therefore, is rooted in the plotting scheme, as in the previous case. 

As an example, we provide in figure [26] images illustrating the ITF quantum phase 
transition: above, T is small and only correlations in the Z-axis are relevant. Below, 
r is large and correlations appear only in the X-axis. The middle panel shows the 
critical case. 

8. Conclusions and further work 

In this work we have described a family of schemes which allow visualization of the 
information contained in quantum many-body wavef unctions, focusing on systems of 
many qubits. The schemes are self-similar by design: addition of new qubits results 
in a higher resolution of the plots. The thermodynamic limit, therefore, corresponds 
to the continuum limit. 

The philosophy behind the schemes is to start out with a region D and divide 
it into several congruent subdomains, all of them similar to D. This subdivision 
procedure can be iterated as many times as needed, producing an exponentially large 
amount of subdomains, each of them characterized by a geometrical index. This index 
can be now associated to an element of the tensor-basis of the Hilbert space, and its 
corresponding wavefunction amplitude goes, through a certain color code, into that 
subdomain. The most simple example is with D a square which splits into four equal 
quadrants, but we can also start with a right triangle, or even with a line segment. 

Physical features of the wavefunctions translate naturally into visual features of 
the plot. For example, within the scheme in section [21 the spin-liquid character of 
the ground state of the Heisenberg model shows itself in a characteristic pattern of 
diagonal lines. This pattern is able to distinguish between open and periodic boundary 
conditions. Other features which show up in the plots is magnetization, criticality, 
invariance under translations or permutation of the qubits, and Marshall's sign rule. 
We have analysed the characteristic features of product states, the ground states of the 
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Figure 26. Frame representation for N = 6 qubits in the ground state of the AF 
ITF Hamiltonian |6] with PBC. The top panel shows the case where F — > 0, and 
correlations are established in the Z-axis. The central panel is critical, F = 1 and 
the lower one shows the F — > oo case, where correlations are strong in the Jf-axis. 



Ising model in a transverse field, the Majumdar- Ghosh Hamiltonian or Dicke states. 
We have also studied spin-1 systems, such as the AKLT state. 

A very relevant physical feature which becomes apparent in the plots is 
entanglement. Factorizability is straightforward to spot: a wavefunction is factorizable 
if all sub-images at a certain division level are equal, modulo normalization. The 
Schmidt rank of a given left-right partition of the system is related to the dimension 
of the subspace spanned by all sub-images within the corresponding subdivision of 
the plot and, so, a crude method to obtain an upper bound is to count the number 
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of different sub-images. The full information about entanglement is contained in the 
matrix that we have termed as cross-correlation, which contains the overlap between 
all subimages at a certain division level. 

In a very different spirit, we have illustrated the frame representations of quantum 
states of many qubits. This approach is related to Wooters' group ideas. In it, the 
expectation values of a selected set of operators are shown in a 2D array, which is 
again displayed in a self-similar manner. 

In this work we have taken the first steps in the exploration of an alternative 
strategy in the study of quantum many-body sytems, which can provide support to 
the corpus of methods in the field. Regarding further work, we would like to stress the 
further exploration of interesting quantum many-body states which we have not done 
here, for example the ground states of fermionic Hamiltonians, the Hubbard model, 
the Mott transition or the BEC-BCS crossover. Understanding the plotting structure 
of matrix product states of low dimension might also result profitable. Moreover, the 
mathematical properties of the mapping itself are worth studying by themselves. 

As a final remark, we would like to announce that source code and further images 
can be found at |http : //qubism. wikidot . com) a webpage dedicated to qubism- related 
resources. 
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